Coherent Population Trapping of Single Spins in Diamond Under Optical Excitation 
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Coherent population trapping is demonstrated in single nitrogen-vacancy centers in diamond 
under optical excitation. For sufficient excitation power, the ffuorescence intensity drops almost to 
the background level when the laser modulation frequency matches the 2.88 GHz splitting of the 
ground states. The results are well described theoretically by a four-level model, allowing the relative 
transition strengths to be determined for individual centers. The results show that all-optical control 
of single spins is possible in diamond. 



Using optical laser fields to manipulate single spins in 
solids is a promising path toward solid-state quantum in- 
formation processing. An important advantage of this 
technique over direct microwave excitation of spin tran- 
sitions is micron-scale spatial resolution, which enables 
selective addressing of individual qubits [1]. Optical spin 
control is also important for interfacing flying and sta- 
tionary qubits as needed for quantum networks [2] and 
repeaters |3[. 

Closely related to this is the effect known as coher- 
ent population trapping (CPT) [4|, observed first in 
gasses [5] and later developed into electromagnetically 
induced transparency (EIT) p . When multiple spin lev- 
els are driven by optical fields to a common excited state 
(a A configuration), a non-absorption resonance can oc- 
cur due to destructive quantum interference between two 
absorption pathways. A dark state forms which is a co- 
herent superposition of two ground states with proba- 
bility amplitudes tunable through the laser amplitudes. 
CPT can be viewed as a steady-state version of optical 
spin control, while time-varying fields allow for dynamic 
control. An important requirement is long-lived ground- 
state spin coherence even under strong optical excitation 
of the material. 

Coherent population trapping and EIT have now been 
obtained in a variety of solids. For example, extremely 
long storage times and room-temperature CPT have been 
achieved in PriYSO @, Q and ruby [9], respectively, but 
the oscillator strengths in these materials seem too small 
for experiments with single impurities. In semiconductor 
systems the oscillator strength can exceed unity, and in 
single charged quantum dots optical pumping [10] and 
initialization of a particular coherent superposition of 
spin states have recently been reported, representing 
a promising step toward all-optical spin control. How- 
ever, in quantum dots as well as shallow donors [12[ the 
spin coherence is thought to be limited by hyperfine in- 



teraction with randomly oriented nuclear spins. Another 
promising system is the nitrogen- vacancy (NV) defect in 
diamond, which has been identified as a promising qubit 
because of its long phase memory [T^j. Here, we show 
that by isolating a single NV center in diamond we can 
obtain a nearly ideal CPT resonance. 

Composed of a substitutional nitrogen next to a car- 
bon vacancy, the NV center can have extremely long- 
lived spin coherence because the diamond lattice is com- 
posed primarily of 12 C, which has zero nuclear spin. The 
NV center has been studied extensively with recent em- 
phasis on potential applications in quantum information 
processing. For the negatively-charged NV center, the 
ground states (Fig. 1) consist of a spin triplet with a 
2.88 GHz splitting between the lower m s = level and 
the upper m s = ±1 levels [14j. These are connected to 
excited states by optical transitions of moderate strength, 
with a total oscillator strength of approximately 0.2 for 
the total vibronic band, or 0.006 for the zero-phonon line 
alone. By exciting these transitions and detecting the 
resulting fluorescence, readout of single spins [15] and 
optically-detected electron spin resonance in a single NV 
center [16| have been demonstrated. This has been ex- 
tended to controlled coupling between the electronic spin 
of a single NV center and a nearby electronic or nuclear 

spin 0, EE E3. 

For NV centers it is not clear in the current literature 
how best to realize a A system due to uncertainty in the 
excited-state structure. Spectral-hole-burning studies on 
large ensembles indicate excited levels which sometimes 



couple to multiple ground sublevels [20|, [2l| , while recent 
work on single centers suggests parallel, spin-conserving 
transitions at zero magnetic field [15]. One approach 
for obtaining a A system is to use a magnetic field to 
mix the ground states [22[. Alternatively, a A system 
can be found even at zero magnetic field due to factors 
such as strain which reduce the symmetry and primar- 
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FIG. 1: Schematic energy level diagram of a strained NV 
center showing resonant optical frequencies coupling ground 
states |1), 1 2) and |3) to excited state |4). are the Rabi 
frequencies, Si are the optical detunings, and a and b are the 
splittings in the lowest three excited states. 



ily modify the excited states [231 ] . The level diagram in 
Fig. 1 is appropriate for the strained case. The excited 
states are split by their orbital components into two spin 
triplets. As discussed in Ref. [24J, under strain the non- 
spin-conserving transitions can be enhanced through the 
spin-orbit interaction. 

We performed measurements on natural type Ha di- 
amond samples obtained from URAL with a (111) ori- 
entated top surface. The samples have low enough NV 
concentrations that single centers can be accessed indi- 
vidually by confocal microscopy [25[. The samples were 
cooled to 2 — 10 K in a liquid- helium cryostat. To increase 
the proportion of NV centers with allowed spin- flip tran- 
sitions, we applied stress mechanically by mounting the 
sample such that it was squeezed by the copper sample 
mount when it contracted upon cooling. A piezoelectric 
transducer allowed for additional adjustments. By this 
method we obtained a splitting in the ensemble zero- 
phonon line of ~ 30 GHz corresponding to a stress of 
~ 30 MPa. Under this condition we could find a spectral 
component for which most of the NV centers had allowed 
spin- flip transitions. However, because of strain and dis- 
order naturally present in many samples, externally ap- 
plied stress is not always necessary to find NV centers 
with allowed spin- flip transitions. The magnetic-field 
measurements described below were performed without 
externally applied stress. 

The main steps in the experiment are shown in Fig. 2. 
An excitation laser, resonant with the 637 nm zero- 
phonon line of the negatively-charged NV center, was 
focused to a diffraction-limited spot within the sam- 
ple. Detection of fluorescence into the phonon sidebands 
from approximately 670 — 750 nm provided a measure 
of excited-state population. Since fluorescence from NV 
centers with allowed spin-flip transitions cannot be ob- 
served with a single excitation frequency due to optical 
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FIG. 2: Experiments at T — 8K with laser modulation side- 
bands approximately 50% as intense as the fundmental com- 
ponent: (a) Fluorescence intensity recorded as the laser fre- 
quency was scanned repeatedly while the modulator was fixed 
at 2.90 GHz. The two bright lines are from a single NV in the 
center of the field, while most of the weak lines are from back- 
ground NVs. (b) Repeated modulator frequency scans with 
laser fixed on m s — transition, (c) Sum of many modulator 
scans, including only those with total fluorescence intensity 
above a threshold. 



pumping, we used an electro-optic modulator to produce 
sidebands at ±2.90 GHz, slightly off resonance from the 
ground-state splitting, to excite all ground states simulta- 
neously. When the laser frequency is scanned, two bright 
fluorescence peaks belonging to the same center are ob- 
served (Fig. 2a). For the left peak, the laser frequency 
excites the m s = ±1 ground states while a modulation 
sideband excites the m s = state. For the right peak, the 
laser excites the m s = state while the other sideband 
excites the m s = ±1 states. Between scans, a pulsed 
repump laser (532 nm) was applied to re-initialize the 
NV center. This was required because the NV center 
was observed to bleach after ~ 10 6 fluorescence cycles of 
resonant excitation. A crucial feature of these diamond 
samples is that single NV centers often exhibit spectral 
diffusion over a frequency range of 100 MHz or less even 
after the repump is applied as in Fig. 2a. 



To observe coherent population trapping, the laser fre- 
quency was fixed on one of the peaks, and the modula- 
tion frequency was scanned across 2.88 GHz with a 10 Hz 
repetition rate (Fig. 2b). Summing over many scans re- 
veals a broad peak centered at 2.88 GHz with a narrow, 
central dip extending almost down to the background 
level, suggesting a large degree of spin coherence on two- 
photon resonance. The blinking between scans is due to 
the stochastic nature of the re-pump process. For the 
data in Fig. 2c, a threshold procedure was used to in- 
clude only those scans (about 50%) for which the NV 
center was in its optically active state. This reduces the 
background level with little effect on the shapes of the 
curves. 

To understand this behavior theoretically, we use a 
four-level model which includes the three ground states 
|1) , 1 2) and |3) representing m s = and two orthogonal 
linear combinations of m s = ±1, and an excited state |4). 
The effective Hamiltonian in a rotating frame under the 
rotating- wave approximation is, 
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where Si and S 2 are the laser frequency detunings from 
the 1-4 and 2-4 transitions, respectively, S 23 is the level 2- 
3 splitting, and fti are the Rabi frequencies proportional 
to the square root of the laser intensities. When spon- 
taneous emission is also taken into consideration, only 
eigenstates with zero probability amplitude in |4) are sta- 
ble (dark states), and all other states decay through flu- 
orescence. Dark states can occur under two conditions: 
f^ll) — ^i|2) is a dark state if Si = S 2 , and ^l) — ^i|3) 
is a dark state if Si = S 2 + £23. For the degenerate case 
S 2 s = we have a two-dimensional dark subspace. To 
obtain a quantitative prediction for the fluorescence in- 
tensity, which is proportional to the excited-state popu- 
lation, we must include spontaneous emission from |4) as 
well as additional decoherence terms. We then solve a 
master equation [26[ for the steady-state density matrix 
p according to, 



dp/dt = -i[H/h, p] + K[p] = . 



(2) 



For TZ[p] we included population decay rates Ti from |4) 
to ground states \i) and decay rates 7^ of the off-diagonal 
density matrix elements pij. 

To test agreement between experiment and theory 
and to obtain information on the individual transition 
strengths, we performed measurements under varied ex- 
citation conditions on another NV center. For the graphs 
on the left side of Fig. 3, the laser was on resonance with 
an m s = transition while a weak modulation sideband 
(about 2% relative power) was scanned across the corre- 
sponding m s = ±1 transitions. In this case the shape of 
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FIG. 3: Measured and fitted fluorescence intensity vs. mod- 
ulation frequency, using the middle excited state of an NV 
with excited-state level spacings (see Fig. 1) a = 0.31 GHz 
and b — 0.94 GHz. (a,c,e): laser fixed on m s = transition; 
(b,d,f): laser fixed on m s = ±1 transition. Excitation powers: 
(a,b): 1.5 ^W; (c,d): 0.5 ^W; (e,f): 0.14/iW. 



the fluorescence dip is determined mainly by the strength 
of the m s = transition. For the graphs on the right, 
the roles of the laser and sideband were reversed. Re- 
sults from several excitation powers are shown along with 
model fits. To reduce the number of free parameters, 
we used the constraints I\ = T = 2n x 13.4 MHz, 

r 3 /r 2 = n 2 3 /nl 7l4 = 724 = 734 = r/2 + 74, and 

712 = 713 = 723 = 71 • For each left/right pair of mea- 
surements, the ratios between and the relevant exci- 
tation powers were held constant. The fits in Fig. 3 used 
Ti/T w 0.8 and 74 « 2tt x 23 MHz and also included 
an adjustable background with linear slope. Under these 
conditions we obtain excellent agreement between the- 
ory and experiment. From the fits we can estimate the 
relative transition strengths for equal excitation power: 
nl/Qj = 0.14 and nj/nj = 0.05. Thus all ground states 
are coupled to the excited state. The m s = ±1 split- 
ting of 5 MHz visible in the graphs on the right could 
be due either to strain or to a background magnetic 
field. We estimate an effective ground-state decoherence 
rate 71 = 2tt x 1.2 MHz, which likely includes optically- 
induced decoherence mechanisms not described explicitly 
in the model. Analysis of the fitted density matrix close 
to two-photon resonance in Fig. 3a shows a nearly equal 
statistical mixture of the two dark states. The ground- 
state coherences pi 2 and pis are approximately 99% and 
90%, respectively, of the maximum possible for such a 
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FIG. 4: Averaged fluorescence intensity vs. modulation fre- 
quency with laser fixed on m s — ±1 transitions for various 
magnetic fields: (a) 0G (b) 4G (c) 10 G (d) 17 G. 

mixture. 

To obtain a definite spin superposition, we can apply 
a weak magnetic field to remove the near-degeneracy of 
the m s = ±1 levels. Fig. 4 shows coherent population 
trapping data on a single NV center for several values of 
magnetic field. The fluorescence dip appears first near 
2.88 GHz modulation frequency for zero magnetic field 
and then splits into two dips following the Zeeman split- 
ting of the m s = ±1 states. Each dip corresponds to for- 
mation of a single dark state. We expect that increasing 
the magnetic field much further would lead to a decrease 
in the fluorescence intensity. Once the m s = ±1 splitting 
is larger than the optical transition linewidth, most of 
the population will collect in the unused ground level. 

These results demonstrate steady-state formation of 
coherent superpositions of spin states of a single NV cen- 
ter with probability amplitudes directly tunable through 
the laser and sideband amplitudes. We expect based on 
these results that dynamic, optical control of the spin 
state will also be possible. Furthermore, this system al- 
lows for a tuning of the relative transition strength of the 
spin-flip versus spin-conserving optical transitions via ex- 
ternal parameters such as strain or electric field. This 
might prove useful as a way to modulate the photon- 
spin interaction. For potential applications, a remaining 
difficulty is how to reduce the spectral instability and 
inhomogeneity of the NV centers. It has recently been 
demonstrated that individual NV centers can be tuned in 
frequency through an applied electric field [251 ]. provid- 
ing a possible solution to this problem. To fully realize 



the potential of diamond NV centers for photonic quan- 
tum information processing it will be necessary to cou- 
ple them efficiently to optical cavities and waveguides. 
This would enhance the modest oscillator strength of the 
zero-phonon line and improve the extraction of emitted 
photons. Then, schemes for efficient interconversion be- 
tween photonic and spin qubits could be realized [27], and 
photons could serve as a communications bus between 
spatially separated spin qubits. Significant progress has 
been made toward fabrication of the necessary structures 
in diamond [28]. 
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